Intermediate DNA at low added salt: DNA bubbles slow the diffusion of short DNA 

fragments 
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We report a study of DNA (150 bp fragments) conformations in very low added salt < 0.05mM, 
across wide DNA concentration range 0.0015 < c < 8 mM (bp). We found an intermediate DNA 
conformation in the region 0.05 < c < 1 mM, by means of fluorescence correlation spectroscopy 
(FCS) and UV-absorbance measurements. FCS detected that in this region DNA has the diffusion 
coefficient, D p reduced below the values for both ssDNA coils and native dsDNA helices of similar 
polymerization degree N. Thus, this DNA population can not be a simple mix of dsDNA and of 
ssDNA which results from DNA melting. Here, melting occurs due to a reduction in screening 
concomitant with DNA concentration being reduced, in already very low salt conditions. The 
intermediate DNA is rationalized through the well known concept of fluctuational openings (DNA 
bubbles) which we postulate to form in AT-rich portions of the sequence, without the strands coming 
apart. Within the bubbles, DNA is locally stretched, while the whole molecule remains rod-like due 
to very low salt environment. Therefore, such intermediate DNA is elongated, in comparison to 
dsDNA, which accounts for its reduced D p . 

PACS numbers: 82.35.Rs 87.15.hp 87.15.-v 



I. INTRODUCTION 

For transcription to proceed in a chromosome, DNA 
has to be separated from the histone proteins and then 
it has to be unwound and the strands separated to allow 
the access of transcription factors (enzymes). This rather 
complicated process has, from a physicist point of view, a 
lot in common with a simple, thermally induced, in vitro 
process of DNA melting/denaturation, where dsDNA he- 
lix completely separates into two ssDNA strands. Sta- 
bility of the double stranded DNA molecule reflects the 
nature of the DNA as a polyelectrolyte 0, H| . When dis- 
solved in polar solvents polyelectrolytes dissociate into a 
polyion and small counterions. Electrostatic interactions 
and entropy of all those charges, as well as those of added 
salt control their phenomenology. Thus, the stability of 
dsDNA depends on added salt, DNA concentrations and 
temperature The DNA melting temperature will 

decrease with a decrease in added salt, and also with a 
decrease in concentration of DNA itself, even more so if 
there is no added salt. The latter is due to a change in 
concentration of counterions. E.g., at high added salt, up 
to physiological levels of the order of 0.1-1 M the DNA 
would melt at high temperatures (80-97°C). Contrary, 
below 1 mM (total counterions and added salt) dsDNA 
may be unstable already at ambient temperature [1, 0] ■ 

These extrinsic parameters are complemented with in- 
trinsic ones: the DNA sequence, conformation and struc- 



ture. Primarily, the sequence of A-T or G-C pairs influ- 
ences the melting. For example, G-C base pairs are con- 
nected with three hydrogen bonds, while the A-T pairs 
with only two, and thus the sequences rich in A-T pairs 
are less thermally stable. The stability of a Watson- 
Crick pair is also determined by its first neighbours, due 
to stacking interactions [8]. This dependence on the lo- 
cal sequence relates to a rather well accepted concept Q 
that the denaturation proceeds from local, cooperative 
openings of several bp long, A-T rich sequences, called 
DNA bubbles, which occur well below melting temper- 
ature [Iol - fl2l |. Close to the melting temperature these 
fluctuations zipper along the strands and the strands sep- 
arate. The bubbles indeed correlate with the actual tran- 
scription initiation sites [l3|. Physically, there is about 
a 10 ksT barrier for opening a bubble and the cost for 
breaking consecutive base-pairs is only about 0.1 fc^T, 
which results in the so-called zippering denaturation pro- 
cess 14j. For the very short fragments, » 10 bp, the 



transition may well be described with a two-state model. 
For a longer DNA molecule, partially open states may be 
conceived, where a part of the sequence is open and part 
still closed. Such molecules, if stable, may be recognized 
to be the intermediate DNA, nor ssDNA, neither dsDNA 
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Models for DNA denaturation start from these general 
concepts. Models may be simple Ising-like, starting from 
the original Poland-Scheraga model [18| , where only the 
energy difference is counted for open and closed base- 
pairs. In more detailed models base interaction depends 
on the distance between the bases (Dauxois-Peyrard- 
Bishop type)[6[. Some models account for the ent ropy 
contribution from the more flexible ssDNA loops (lj|l9|. 
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and find a continuous (II order) or discontinuous (I or- 
der) transition [2(| [HJ . The latter result was crucially 
based on taking into account the self-avoiding interac- 
tions between the various parts of the chain. 

For the experimentalist, it is not straightforward to 
demonstrate the existence of the DNA bubbles. Cer- 
tainly, experiments evaluating the basepair opening prob- 
ability rates and bubble lifetimes, in any specific set of 
conditions for their existence, would contribute to the 
modeling effort and our general knowledge of DNA sta- 
bility and dynamics. The information on base-pair open- 
ing was initially provided by NMR study [ll|, [22| of the 
imino-proton exchange, i.e. study of protons, formerly 
participating in hydrogen bonds, that bases exchange 
with D 2 water when a base pair opens. This NMR 
study actually points that the bubble opening is a very 
local phenomenon, occuring on the level of a single base- 
pair. Opening probabilities (dissociation constants) at 
25°C are actually quite low: 1CT 5 - 1CT 6 for an AT pair 
and 10~ 6 — 10~ 7 for a GC pair. The lifetime of a base- 
pair is of the order of 10 ms and the rate for closing 
the bubble is 10 8 s -1 . Only recently these NMR stud- 
ies have been complemented by other studies of bubble 
dynamics. An FCS study by Libchaber et al. [H| has 
found that fluctuations occur with closing rates 3-4 or- 
ders of magnitude higher than estimated by NMR, i.e. 
10 4 - 10 5 s _1 . The size of the bubbles was 2 to 10 bp 
on a fragment only 18bp long. These authors found no 
contradiction with NMR rate estimates, concluding that 
NMR registers much smaller conformational change with 
respectively higher fluctuation rates. That is, they con- 
sidered that the openings detected by FCS are more bi- 
ologically relevant as they relate to larger, nm-scale sep- 
arations. 

Peyrard et al. [It], [H| were able (by two-photon excita- 
tion of guanines) to detect when GC pairs were not paired 
in short ss 50bp DNA while temperature was raised and 
DNA melted. Thus they obtained localized informa- 
tion on bubble opening and provided sequence dependent 
DNA melting curves. Importantly, they noted that short 
AT rich sequences (TATA boxes) would fluctuate and 
disturb both neighbouring and distant GC pairs well be- 
low the melting temperature. The fact that the opening 
fluctuations can have a non-local (non-single basepair) 
effect is related to nonlinear dynamics of DNA [9] , but it 
was emphasized that this effect becomes negligible as it 
averages out for kbp or longer DNA sequences. 

From the above it may be recognized that intermedi- 
ate DNA states should exist below the melting transition, 
however only one group of authors has insofar managed 
to quantify the fraction of intermediate DNA in the pop- 
ulation of DNA molecules. Montrichok et al. [25| have 
designed DNA constructs, ones with the TATA box at 
one end, and the others with TATA box ("soft" region) 
in the middle and GC clamps at both ends. When melt- 
ing was stopped before complete denaturation, and the 
solution cooled quickly (quenched), strands which were 
completely separated were not able to find their comple- 



ments and formed hairpins instead. Their fraction was 
quantified by simple gel electrophoresis. Fraction of the 
open molecules was compared (at equal temperatures) to 
the fraction of open basepairs in the solution, quantified 
by UV-absorption measurements. If the latter fraction 
was higher then there must have been some partly de- 
natured, intermediate molecules before the quench. In- 
termediate DNA states identified from these experiments 
were succesfully modeled by Monte Carlo simulations us- 
ing the Dauxois-Peyrard-Bishop model [l6| . 

One issue that arises is whether the local fluctua- 
tions influence mechanical properties of the whole DNA 
molecule. Early experiment, where self-diffusion coef- 
ficient of 17 kbp DNA molecules was probed in tem- 
perature by dynamic light scattering found no effect of 
fluctuations on DNA rigidity (persistence length) [26[. 
Only recently it has been noted that there is the effect 
of fluctuations, but only on very short (10 -90 bp) DNA 
fragments, which could behave as having a persistence 
length as low as 10-20 nm p7, 28], much shorter than the 
agreed value L p — 50 nm [29j |. Also, Tomic et al. have 
noted that persistence length of DNA might be reduced 
in some conditions [30( ] - They worked with mononucleo- 
somal w lbObp DNA fragments, which have the contour 
length L c = 50 nm similar to the persistence length, so 
they are practically rodlike. Their dielectric spectroscopy 
technique, which provides characteristic length scales of 
a polyelectrolyte system, complementary to e.g. data by 
small angle X-ray scattering [31] . indeed found a 50 nm 
length scale for fragments in pure water. However for 
fragments in 1 mM NaCl the characteristic length scale 
was 20-30 nm, lower both than L c or L p . This would 
indicate that added salt, which does reduce the electro- 
static contribution to the persistence length [32j, has re- 
vealed L p lower than the agreed value. We may presume 
that the reduced L p is the consequence of partial DNA 
denaturation in thesee conditions of low salt. 

In this paper, we propose that intermediate DNA is the 
major constituent in mononucleosomal DNA (denoted as 
DNA146) conformations population in the concentration 
range 0.05 < c < 1 mM (bp) at very low salt conditions, 
Csait < 0.05mM. We infer the presence of intermediate 
DNA from the DNA146 polyion self-diffusion coefficient 
D p , obtained by the fluorescence correlation spectroscopy 
measurements (FCS, method described in Sec. II). In 
Sec. Ill we present the D p value for 0.05 < c < 1 mM 
range which is below the values for both ssDNA and na- 
tive dsDNA. This means that DNA conformation pop- 
ulation there can not be a simple mix of native dsDNA 
helix and separated ssDNA strands. In Sec. IV we argue 
that an intermediate DNA conformation must be present, 
hydrodynamically larger than both ssDNA or dsDNA of 
similar polymerization degree. We suggest that under 
these conditions the basepair opening and closing rates 
have changed in such a manner to induce a variation of 
the diffusion properties of the whole molecule. In this 
sense, the intermediate DNA may be considered to be a 
stable form. Eventually, we analyze D p data and extract 
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the fraction of completely open DNA molecules, ssDNA, 
against partly open, intermediate DNAs and native ds- 
DNA. We also compare these fractions with measured 
fraction of open DNA base-pairs, following the procedure 
of Montrichok et al. [25| • 

II. MATERIALS & METHODS 
A. Samples 

We studied solutions of nucleosomal dsDNA fragments 
about 150 bp and 50 nm long. This DNA is denoted 
as DNA146, since there are 146 bp of DNA wrapped 
around a histone octamer in a nucleosome core particle. 
DNA146 pellets from the same stock as in [33[ was used. 
The monodispersity of fragments was checked by gel elec- 
trophoresis. We also checked whether the pellets contain 
any added salt. DNA concentrations will be expressed 
as molar concentrations of basepairs in the remainder of 
this paper. 

Several solution sets with different DNA146 and added 
salt concentrations were prepared, which we describe be- 
low: 

set I Solutions with varying DNA concentration, in very 
low added salt: A set of 40 solutions covering the 
range of DNA concentrations c = 0.001 — 8mM was 
prepared by dilution with pure water of aliquots 
from the mother solution to reach the desired con- 
centration. The mother solution was obtained by 
dissolving 10 mg of dry DNA146 in 0.55 mL pure 
water - the final concentration was c = 27 mM. The 
UV-absorbance at 260 nm was measured, at 25°C 
for all the solutions in the set. Only then, for FCS 
measurements, the DNA146 solutions were fluores- 
cently labeled by addition of synthetic HObp DNA 
from a 0.5 mM stock kept in 10 mM TE buffer. This 
HObp dsDNA will be further denoted DNA110*. 
The asterisk * denotes the fluorescent labelin g by 
Cy5 fluorophore at one end of DNA molecule [33| ■ 
DNA110* was added to a concentration 1.5-2 fjM 
to achieve 20 nM Cy5 label concentration. The 
buffer of DNA110* stock gets diluted after addi- 
tion into pure water DNA146 samples, thus intro- 
ducing a minimal added ionic strength, not above 
0.05mM. In this manner we regarded these samples 
as very low salt solutions, c sa it < 0.05mM, and la- 
beled them appropriately. FCS results for this set 
are presented also in |33j . 

set II Solutions with varying DNA concentration, in 10 
mM TE buffer. From c = 27 mM pure water 
mother solution (see above) a 0.1 mL sample was 
taken and diluted with 0.25 mL of 35 mM TE 
buffer. We used 10:1 TrisHClEDTA buffer set to 
pH7.5. Thus we obtained DNA146 starting so- 
lution with c = 7.7 mM in 10 mM TE buffer . 
For FCS measurements DNA110* was added to 



this solution, to a concentration 1.5-2 /zM. After 
FCS measurement the sample was further diluted 
with appropriate volume of 10 mM TE buffer to 
produce a next lower DNA146 concentration sam- 
ple. The procedure was repeated and 15 different 
DNA concentrations were measured in the range 
c = 0.03 — 8mM. On each dilution appropriate vol- 
umes of DNA110* stock were being added to main- 
tain the 20 nm Cy5 level. 

set III Solutions with fixed DNA concentration, in varying 
added salt: DNA146 mother solution, c = 7.5 mM 
was prepared by dissolving a pellet in 10 mM NaCl. 
By dilution with appropriate volumes of pure water 
and ImM NaCl we prepared a solution set where 
DNA concentration was always c = 0.75 mM, while 
NaCl concentrations were 0.2, 0.5, 1, and 2 mM. 



B. UV absorption measurements 

DNA146 sample set I (varying DNA, very low salt 
< 0.05 mM) has been tested without dilution for UV- 
absorbance at 260 nm with Nanodrop (ThermoScien- 
tific). This instrument shows a linear response in ab- 
sorbance for a very broad DNA concentration range, 
c = 0.003 — 15 mM, with the sample volume of only 
1-2 /iL. The Nanodrop instrument was critical to get the 
actual UV absorbance of DNA samples without dilution. 

For DNA melting experiments on set III (fixed DNA, 
varying salt) in 0-90 °C range we have used Agilent - 
HP 8452 Spectrophotometer with a Peltier based tem- 
perature control. From the temperature dependent UV- 
absorbance at 260 nm (hyperchromicity) curves we ex- 
tracted the melting transition temperature T m and tran- 
sition width dependence on total ion concentration. 



C. Fluorescence correlation spectroscopy 

We have used Zeiss ConfoCor II FCS instrument. 
Focal volume was defined by a Zeiss Plan-NeoFluar 
100x/NA1.3 water immersion objective, epi-illumination 
was by He-Ne 632.8 nm 5mW laser, for excitation of Cy5 
fluorophore. Measurements were performed at 25°C, the 
ambient temperature of the temperature stabilized clean- 
room. 

Fluorescence correlation spectroscopy is used to mea- 
sure the diffusion coefficient of the fluorescently labeled 
molecules. Number fluctuations of the molecules enter- 
ing and leaving the focal volume of the instrument are 
registered as fluorescence variation. An autocorrelation 
function is calculated for the fluorescence intensity trace. 
This function decays exponentially with autocorrelation 
time r c : 
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G ^ = F' TT^ / ^==(1+ T ^ f «cp(--)) 

JVf 1 + "T Ml | ( w 0)2T e \ 1 — J TT 

(1) 

iV/ is the average number of fluorescent molecules in the 
focal volume, zq/wq is the focal volume structure pa- 
rameter and T, average fraction of fluorophores in the 
triplet state (thus, non- fluorescing) . The lifetime tt of 
the triplet state is taken into account when fitting. The 
characteristic decay time r is the diffusion time that the 
fluorescent molecule takes to traverse the focal volume. 
The details of the procedure we used to extract r may 
be found in [331 ] - 

In FigfTJ besides the autocorrelation function experi- 
mental curve and fit for Cy5, we present the G(t c ) curves 
recorded for 2 of DNA110* found in either 0.2 mM 
DNA146 in very low salt c salt < 0.05mM, (set I) or 0.2 
mM DNA146 in 10 mM TE buffer (set II). The experi- 
mental curves have been normalized, and thus they ap- 
pear to overlap due to a rather small difference in the 
diffusion times. Nevertheless, the diffusion times r that 
we extracted are different, with values of 390 and 430 
/is, respectively. To demonstrate that the fits reliably 
distinguish between the two experimental data sets for 
DNA110*, the area of the main panel, denoted with the 
rectangle, is shown enlarged in the inset [34j]. The ar- 
row in the inset denotes the 40 /its difference between the 
respective diffusion times obtained by fits. 

The diffusion time is inversely proportional to the self- 
diffusion coefficient D p of the molecule. The relationship 
between r and D p may be obtained from the measure- 
ment of Tc y 5 for Cy5 fluorophore, whose diffusion coeffi- 
cient is known, D Cy5 = 3.16 • 10" 10 m 2 /s HU: 



D p = D c 



y5 " 



T C y5 



(2) 



The diffusion times D p that we can obtain from r are for 
the fluorescently labeled DNA110* molecules at a rela- 
tively low concentration, diffusing in solutions of varying, 
but mostly higher DNA146 concentration. That is, we 
obtain D^^, diffusion coefficient for DNA110* which, 
however, depends on DNA146 concentration c. However, 
we assume that it is possible to rescale ^iio*( c ) to get 
the values for DNA146 D^(c). We remind that dsDNA 
persistence length L p = 50 nm [29| is comparable to con- 
tour lengths of both DNA110* and DNA 146, 38 and 50 
nm, respectively. Therefore, we expect that they assume 
an extended, rodlike configuration. According to Tirado 
et al. [36[ the translational diffusion coefficient calculated 
for a rodlike macromolecule is given by 



D 



th 



kT \n{L c /d) + 0.312 
3irr) L c 



(3) 



Here L c = Nb is contour length, d is polyion diameter, 
T) = 8.9- 10~ 4 Pas is viscosity of water (T = 298 K). With 
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FIG. 1: Experimental autocorrelation functions G(r c ) (de- 
noted DNA110*) recorded for a sample from set I (2 /xM 
DNA110* in 0.2 mM DNA146 in very low salt < 0.05mM) 
or a sample from set II (2 fiM DNA110* in 0.2 mM DNA146 
in 10 mM TE buffer). For comparison, the autocorrelation 
function (denoted Cy5) is shown for 20 nM Cy5 fluorophore 
in pure water, without any DNA. The inset zooms the corre- 
lation time range 250 — 550/is, denoted by a rectangle in the 
main panel. The experimental values of G(t c ) are shown with 
symbols and the respective fits with lines. 



b = 0.34 nm and d = 2.6 nm, the diffusion coefficient for 
HObp DNA is Dff „ = 3.98 • 10" 11 m 2 /s and for 146bp 
DNA is Z^V = 3.27 ■ 10" 11 m 2 /s. Stellwagen et al.^ 
have reviewed the literature and shown that the expres- 
sion by Tirado et al. is well applicable to experimental 
data obtained for DNA molecules in size from 10 to 1000 
basepairs. 
Thus, 



D?£(c) 



pith 

-^110* 



(4) 



In this manner, fluorescence correlation spectroscopy was 
used to obtain the self-diffusion coefficient of DNA146, 
DlH(c) for varying DNA (DNA110* and DNA146 com- 
bined) concentrations c = 0.0015 — 8 mM. Most impor- 
tantly, we were able to distinguish and compare the val- 
ues obtained for samples in very low salt (set I) and sam- 
ples in 10 mM TE buffer (set II). 



III. RESULTS 

A. DNA melting in very low salt conditions 

Classically, UV absorbance measurements at 260 nm 
are used to study the DNA melting, i.e. helix-coil tran- 
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FIG. 2: Melting curves for DNA146, c = 0.75 mM (bp) 
with varying added salt 0.2, 0.5, 1, and 2 mM NaCl (set III), 
denoted as a, b, c, d, respectively. Melting appears as an 
increase in cuv/c ratio. The arrows denote melting tempera- 
tures T m defined at half height of the maximum cuv/c rise. 



sition 0, 0] • Since the extinction coefficient for the ds- 
DNA, 1 O.D.=50 /ig/ mL is about 40% higher than for 
ssDNA, 1 O.D.=35 /zg/mL, concomitantly with temper- 
ature induced DNA denaturation, occurs a 40% rise in 
absorbance (hyperchromicity effect). 

In FigJ21 we present thermal denaturation curves for 
DNA146 samples at a fixed DNA concentration c = 
0.75 mM, in varying added salt, 0.2-2 mM NaCl. The 
measured UV absorbance rise is presented as the rise 
of the ratio between UV-spectrophotometrically deter- 
mined concentration cjjv and the gravimetrically deter- 
mined (during sample preparation) concentration. The 
absorbance was converted into cjjv concentrations us- 
ing the extinction coefficient for the dsDNA. cjjy /c ratio 
close to 1.4 thus corresponds to fully denatured DNA and 
cuv/ c = 1.0 corresponds to the native dsDNA. Melting 
temperatures T m are therefore defined at half height of 
the maximum cjjv/c rise, where half of the basepairs 
should be open. For DNA146 in 0.2 mM NaCl solu- 
tion we note that already at ambient temperature UV- 
absorbance, (cjjv/c > 1) is enhanced. Therefore we pre- 
sume that this sample was already partly denatured at 
ambient conditions, due to the low salt content. 

In Fig. [3] we present the ratio cjjv / c, measured at 
25°C for DNA146 dissolved in very low added salt 
Csait < 0.05 mM (set I). At higher concentrations, above 
c = 1 mM cjjv and c seem to comply indicating that 
the DNA remains native. At the lowest concentrations 
cjjv I c increases to about 1.4, which easily relates to 40% 
hyperchromicity of denatured DNA. Therefore, in Fig. [3] 
dashed lines at cjjv /era 1.4 and 1 are respectively labeled 
ssDNA and dsDNA. The scatter in the data is mostly 
due to errors in DNA sample preparation, and to some 
extent due to the very minute sample droplets applied 
(a cuvette free spectrophotometer is used), which may 
evaporate and change concentration rather quickly. 
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FIG. 3: The cuv/c ratio for DNA146 in very low salt c sa it < 
0.05 mM (set I). The dashed lines at cuv/c ~ 1.4 and 1 denote 
fully denatured ssDNA and native dsDNA, respectively. The 
arrows and letter e denote the concentration range where the 
ratio cuv/c is just below 1.2, halfway between ssDNA and 
dsDNA levels. 



Further, we note DNA146 samples from the interme- 
diate concentration range c = 0.4 — 0.75 mM that show 
cjjv I c ra 1.2, corresponding to about 20% hyperchromic- 
ity (c/. Figf2]). These samples are thus halfway through 
the melting transition, at ambient temperature. Thus, 
for DNA146 in the c = 0.4 — 0.75 mM range we estimate 
the temperature where we performed measurements to 
be the melting temperature, T m — 25°C. 

In FigH] we show DNA146 melting temperature T m 
vs. total ionic concentration of the solutions. We show 
the values obtained from melting curves (sample set III, 
FigUl points a, b, c, d) alongside the just above defined 
T m = 25°C for DNA146 c = 0.4 - 0.75 mM concentra- 
tion range (sample set I, Fig 121 points e). Total ionic 
concentration ctot is a sum of all the ions present in the 
solution: counterions c, — 2c and added salt anions and 
cations 2c NaC i = c Na + + c ci - , thus c to t = <k + 2c Na ci- 

As predicted by Manning 38| we obtain the linear 
slope dT m /d\ogc to t — 28 ± 3°C per decade, which com- 
pares very well with the slope 27°C per decade found 
by Record [4] for T4 (170 kbp) and T7 (40 kbp) bacte- 
riophage DNA solutions without added salt. The slope 
reflects the counterion condensation and should be inde- 
pendent of the chain size and thus invariant from the type 
of DNA used in the experiment, as, indeed, we observed. 
While the slope is similar, we note that the transition 
temperatures for DNA146 are about 10°C lower than for 
long T4 and T7 DNAs, for equal Ctot- We also note that 
the melting transitions for DNA146 are much broader, 
ra20-30°C, than for T4 and T7 DNAs, where they are 
only 4-9°C. That the melting temperature decreases and 
the transition width increases for shorter DNA chain size 
is to be expected [H, Ho| ■ 

The above determined ability of DNA146 to keep the 
native conformation at ambient temperature, in the very 
low salt (< 0.05 mM) environment, but only above a cer- 
tain DNA concentration c = 1 mM. reflects the findings 
by Record [4| and by Tomic et al. 4l|. In the following 
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FIG. 4: DNA146 melting temperatures T m , as a function of 
total ionic concentration a + 2cNaCi of the solution. Coun- 
terion concentration Ci is denned as twice the DNA base- 
pair concentration c. Circles denote T m data obtained for 
c = 0.75 mM DNA146 in varying salt (points a, b, c,d, set 
III) solutions. For DNA146 c = 0.4 - 0.75 mM samples in 
very low salt (points e, set I) we have estimated T m = 25°C 
(triangles). 

we identify the DNA conformations present in the very 
low salt solutions, with the region below c = 1 mM being 
the most intriguing. 



B. DNA146 diffusion coefficient 

In Fig. [5]we present -D^g, the diffusion coefficient of 
DNA146 polyion as a function of DNA concentration c 
in the range 0.0015-8 mM. The coefficient values were 
derived from diffusion times measured by FCS, (see Ma- 
terials & Methods). Both D%$&c) for DNA146 in very low 
salt < 0.05 mM (sample set I) and for DNA146 in 10 mM 
TE buffer (sample set II) are shown. The concentration 
range of the FCS study encompasses the range studied 
by UV-spectroscopy. According to behavior of D^(c) 
the studied concentration range has been divided into re- 
gions A,B,C,D with respective delimiting concentrations 
c AB = 1 mM, c B c = 0.05 mM and c CD = 0.006 mM. 

In Fig. [5] we also show the diffusion coefficients D\\ 6 
and -Df|g, extrapolated for a DNA of a polymerization 
degree N = 146. These are extrapolated (according to 
Eqs l2l4p from the diffusion times measured for 1.5 /jM 
(basepair) of fluorescently labeled N = 110 ssDNA and 
native dsDNA (DNA110*), respectively. First, to get 
-Df! 6 we used a sample of fluorescently labeled, single- 
stranded, 110 nt long oligomer, in pure water solution. 
The diffusion time was extrapolated according to Eqs l2l4l 
and wc got Df| 6 = 3.95 • 10~ n m 2 /s. Then, to get D% 
we used DNA110* added to 10 mM TE buffer, with- 
out any DNA146. Since in buffer, DNA110* is certain 
to keep the native, double helix conformation. Thus, 
-Dffe = 3.2 • 10 -11 m 2 /s is considered as the represen- 
tative experimental value of the diffusion coefficient for 
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FIG. 5: Diffusion coefficient D^i for DNA146 in very low 
salt < 0.05 mM (inverted triangles, sample set I) and in 10 
mM TE buffer (triangles, sample set II). D^j- varies with 
the DNA concentration c (basepair) across regions denoted 
A to D. Full line denotes the fit to EqO which describes the 
concentration dependence of diffusion coefficient for rod-like 
particles. Lower dashed line shows the theoretical diffusion 
coefficient value D{\ & = 3.27 • 10 -11 m 2 /s for native DNA146 
helices in dilute solution. It also points to the value Dfl 6 ex- 
trapolated from the diffusion time measured for DNA110* in 
lOmM TE buffer, without DNA146 (see text). Upper dashed 
line points to the diffusion coefficient value (inverted black 
triangle)DiI 6 = 3.95 ■ 10 -11 m 2 /s expected (see text) for 146 
nt ssDNA strands in pure water. 



native dsDNA, N = 146. 

First we note that, in region A (c = 1 — 8 mM), the 
values DlH(c) for DNA146 both in 10 mM TE buffer 
and in very low salt < 0.05 mM, are equal within the 
experimental error. This result shows that DNA is of na- 
tive conformation in both environments. The variation 
of is due to the 50 nm DNA146 fragments being in 
the semidilute regime above 1-2 mM [3fJ|. That is, the 
polyion diffusion slows down, due to an interchain repul- 
sive interaction that increases with reduction of average 
distance between polyions, i.e. with the increase in con- 
centration of polyions. Self-diffusion coefficient, D p (c) 
has been previously shown to decrease linearly with the 
volume fraction $ of rod-like particles, for e.g. a mineral, 
boehmite, L w 300 nm, 42] and 20 bp dsDNA L w 7 nm, 

D p (c) = D° p (l - a*) (5) 

Here, D p is the self-diffusion coefficient measured in di- 
lute solution. For DNA, molar concentration c is related 
to volume fraction as $ w 0.001c. Proportionality factor 
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a depends on the aspect ratio L/d of the particle [42j . 
For DNA146, p « 20 - 25, leading to a « 15. Fit of EqE] 
to our DNA146 in TE buffer data gives a » 30 (full line 
in Fig. [SJ, in a reasonable agreement with the theoreti- 
cal value. Measurements at higher concentrations would 
be necessary to elucidate significance of this discrepancy, 
however this is out of scope here. 

Below cab = 1 mM (basepair), DI^q(c), when mea- 
sured in 10 mM TE is practically concentration indepen- 
dent across almost two decades in concentration (regions 
B,C). These values are very close both to just above 
defined experimental value for 146 bp dsDNA Dfl 6 — 
3.2-10" 11 m 2 /s, and to the theoretical D\\ & = 3.27-10- 11 
m 2 /s value, calculated according to Eqj3]for a rodlike 146 
bp dsDNA, 50nm long with d = 2.6 nm diameter. That 
is, DNA146 in lOmM TE buffer keeps the native state. 
Also, since this is dilute regime, diffusion coefficient is 
constant. 

On the contrary, (c) in very low salt, deviates from 
the native DNA146 value, below cab = 1 mM (basepair). 
Specifically, in the intermediate concentration region B 
(c = 0.05 - 1 mM) = 2.9 • 10" u m 2 /s is rather 

constant, but below (beyond the experimental error) the 
native DNA146 value. This reduction in D^ic) means 
that in region B, DNA polyions diffuse more slowly, than 
DNA polyions which keep the native conformation. 

Below cbc — 0.05 mM (basepair), in region C, D e {^ 
measured in very low salt starts to rise continuously and 
overpasses the native dsDNAl46 value. Eventually, be- 
low ccd = 0.006 mM, in region D, finally rises 
sharply to come very close to -Dff 6 = 3.95 • 10~ n m 2 /s, 
the value that we ascribe to 146 nt ssDNA. Here we 
note that Dfl 6 value is higher than the theoretical value, 
3 • 10~ n m 2 /s for a rodlike ssDNA that may be calcu- 
lated according to EqO using b — 4.3 nm and d = 1.2 
nm. Further, D\\ & value is lower than the experimental 
value 5 • 10 -11 m 2 /s that may be extrapolated for 146 nt 
ssDNA in the high salt (10 mM or higher) j44|. In these 
high salt conditions ssDNA has a persistence length of 
the order of 1 nm and assumes a conformation of a com- 
pact, random walk coil [H, |H|. E.g., 146 nt ssDNA, 
which has a contour length of about 60 nm, would form 
a coil with a hydrodynamic radius of 3 nm. Considering 
the hereabove presented range of diffusion coefficients for 
ssDNA, we judge that ssDNA in our, very low salt con- 
ditions neither assumes the form of a compact coil, nor 
a rigid rod, but a rather extended coil. 



IV. DISCUSSION 

A. Fraction of open basepairs 

As we have noted in previous Section, FCS measure- 
ments show that DNA146 in very low salt < 0.05 mM 
(sample set I) assumes native, double-helix (and rodlike) 
conformation, but only in region A, above DNA concen- 
tration cab = 1 mM. Also, we found that in region D, be- 



low ccd = 0.006 mM DNA appears to be single-stranded, 
forming a rather extended coil. In this manner, we rec- 
ognize the concentration of DNA itself as a parameter to 
induce DNA melting in very low salt environment. 

Therefore we propose that our UV-absorbance data, 
that is, cjjv/c ratio presented in FigSl may be regarded 
as a DNA melting curve, recorded, however vs. DNA con- 
centration as the parameter which induces melting. We 
remind that UV-absorbance only provides the amount 
of basepairs that have opened in a given DNA sample. 
Therefore, analysis of cuv/c data should provide the 
fraction f{c) of open basepairs in the ensemble of DNA 
molecules in samples from set I. 

We also use the Manning free counterion fraction 9(c) 
data that we presented in Ref. [33j , in support of cuv / c 
data. Both data sets have been obtained for sample set I. 
The reported 9(c) variation occurs within the range de- 
fined by the theoretical 9 values for dsDNA and ssDNA. 
The relative variation that we obtain is also similar to 
the relative difference in 9 values for dsDNA and ssDNA 
observed experimentally by other authors [7]. Thus, we 
consider that 9(c) curve represents DNA melting in low 
salt, and measures the fraction of open basepairs, as well 
as cuv I c data. 

The open basepair fraction is derived in a following 
manner 

/(c) = c) z A : m \ (6) 

Here A m i n and A max are the minimum and maximum 
value for the parameter being analyzed, either A(c) = 
cuv /c or A(c) = 9(c). 

Fig. [6] shows a very good correlation between the 
open basepair fraction calculated from UV-absorbance 
f(cuv/ c ) an d from free counterion fraction f(9(c)), 
across the studied DNA concentration range c = 0.015 — 
8 mM. We judge that the f(c) that we obtained is, indeed, 
a very god measure of fraction of open basepairs present 
in very low salt solutions of varying DNA146 concentra- 
tions. Both cuv /c and 9(c) data appear to be valid for 
characterization of the melting transition. 

B. Intermediate DNA 

The most important result of this study is that in 
region B the diffusion coefficient D^q measured for 
DNA146 in very low salt is diminished compared to val- 
ues both for ssDNA and native dsDNA (DNA in 10 mM 
TE buffer). In this same region a fraction between 10 
and 90% basepairs are open (cf. Fig|6]and FigJSJ. 

Careful consideration of the above leads to a conclusion 
that the DNA population in region B is NOT a simple 
mix of native dsDNA helix and separated ssDNA strands. 
If it were a mix then the measured diffusion time and the 
corresponding diffusion coefficient should be of some in- 
termediate value [451 ] . between the experimental values 
for native dsDNA and ssDNA, Df| 6 = 3.2 ■ 10~ n m 2 /s 
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FIG. 6: The fraction / of open basepairs vs. DNA146 base- 
pair concentration. The circles are derived from cuv/c ratio 
(see also FigEJ) . For comparison we show data (diamonds) de- 
rived from free counterion parameter 8(c) taken from Ref. [33| . 
Fraction / » above 1 mM (region A) relates to dsDNA, 
while / approaching the value of 1 below 0.05 mM (region C) 
means that only ssDNA remains in the solution. 



and -Df| 6 = 3.95- 10 -11 m 2 /s, respectively. However, the 
actual value is only 2.9 ■ 10 -11 m 2 /s! The DNA molecules 
there diffuse more slowly than either a rigid-rod dsDNA 
or an ssDNA strand. In other words, the DNA146 con- 
formation present in this region is distinguished by being 
hydrodynamically larger entity than native DNA146 he- 
lix. As this is a dilute solution regime, the interactions 
between the molecules are not a plausible explanation, 
as they are for the decrease of (c) above 1-2 mM, in 
the semidilute regime (region A, FigJS]). 

Wilk et at have also found a diffusion coefficient to be 
reduced for DNA fragments (20 bp) in low salt, 0.01 mM 
NaCl in comparison to high salt, 0.2 M NaCl solutions 
and for a rather similar DNA concentration range (43[. 
However, they have neglected in their analysis the fact 
that DNA denatures in the conditions of very low salt, 
which we have carefully analyzed in the previous sub- 
section. To account for the reduction in diffusion ceoffi- 
cient, they have proposed that the hydrodynamic volume 
of the DNA molecule is enlarged, by the amount of the 
Debye screening length that is, the contour length 
L c becomes L c + k _1 , and the diameter becomes d + n^ 1 . 
We consider this scenario to be inappropriate. Already 
for DNA146 in 10 mM TE buffer (and a rather small 
= 1 nm), the value of the diffusion coefficient D['l 6 
thus calculated (using Eq(3]) would deviate for 10% from 
the experimental result we got. We remind, we get the 
discrepancy of only 1-2% between our experimental Df% 6 
result and theory, if we do not use a correction like Wilk 
et ai. 

Theoretically, the diffusion coefficient for rodlike 
molecules (see EqJ5]) is linearly dependent on the con- 
tour length of the molecule, and only logarithmically on 
its diameter. Thus, the most plausible mechanism for 
D^2g(c) becoming reduced in region B would be if DNA 
contour length increased about 10%, leading to a 10% de- 



crease in diffusion coefficient. Now, we note that dsDNA 
monomer size is 0.34 nm while for ssDNA it is 0.43 nm, 
20% longer [44|. When a sequence of several basepairs 
in a dsDNA molecule breaks hydrogen bonds and un- 
stacks (DNA bubble), then each strand locally assumes 
an ssDNA coil conformation. Consequently the molecule 
would elongate in the region of the bubble. The molecules 
being elongated by 10% - half the maximum 20% - should 
then have half of the basepairs open. Simply, half of the 
sequence would be DNA bubbles and half native helix - 
an intermediate DNA would form. We assume that the 
fraying, separating ends of a dsDNA molecule would have 
a similar effect as DNA bubbles, on slowing the interme- 
diate DNA diffusion. 

Remaining issue would be whether these elongated, in- 
termediate DNA molecules really maintain rodlike shape, 
otherwise application of EqJ3] would not be justified. We 
remind that, in these very low salt conditions, the charges 
along the DNA backbones repel strongly due to the lack 
of Debye screening and it is very plausible that this al- 
lows the intermediate DNA molecule to remain extended, 
rodlike in shape, despite having defects of lower intrinsic 
rigidity |27H30T l4ll . l44j . We provide an illustration above 
the data panel of FigEJ 



C. Fractions of DNA conformations 

For the analysis of fractions of conformations in the 
population of DNA146 molecules, we will postulate that 
intermediate DNA146 is of a rodlike form and is char- 
acterized with about 50% basepairs open, and that the 
only other conformations that may appear in here ana- 
lyzed solutions are fully separated ssDNA coils and rod- 
like native dsDNA helices. In other words, we do not ex- 
pect DNA intermediate states with varying proportions 
of open base-pairs to coexist along the melting transition. 
Accordingly, we will analyze the diffusion coefficient data, 
and find how the fractions of dsDNA helices, intermediate 
DNA and strands of ssDNA vary with the DNA concen- 
tration, in very low salt solution, Figj7l Eventually, we 
will provide arguments to support this proposition. 

First, we compare diffusion coefficient data obtained 
for DNA146 in very low salt (set I) and in buffer (set II) 
m regions A and B. We denote these D%$(c) and Df s i6 (c), 
respectively. Specifically, Dfl 6 (c), being obtained for 
DNA in buffer is thus the value for native helices, where 
the concentration dependence is due to samples above 
c = 1 mM being in the semi-dilute regime. Variation in 
£>i46 (c) is also due to this influence, but, as DNA concen- 
tration is lowered below about 1-2 mM, £)^g (c) deviates 
from Dfl 6 (c), showing that intermediate DNA starts to 
appear in the molecular population. The difference be- 
tween £)i4g (c) and Dfl 6 (c) stabilizes in the lower part of 
region B, c = 0.05 — 0.3 mM, indicating that only inter- 
mediate DNA remained in the population. We remind 
that the diffusion coefficient value D"£ 6 — 2.9 • 10 
m 2 /s in this region is lower than values either for ds- 
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FIG. 7: Fraction of different conformations in the popula- 
tion of DNA146 molecules varies with DNA concentration (in 
basepair) in the very low salt (< 0.05mM) solutions. Closed 
symbols denote the intermediate DNA fraction im(c). Open 
symbols denote the other conformation present in the popula- 
tion in the respective region (1 — im(c)). Circles and squares 
(both open and closed) denote fractions obtained according to 
Eq[7]&[8l respectively. In region A, dsDNA dominates, start- 
ing to be replaced by intermediate DNA towards region B. 
Across most of the region B, intermediate DNA is dominant 
in the population. In region C intermediate DNA is replaced 
by ssDNA. Eventually, in region D, ssDNA dominates. Above 
the data-panel DNA conformations are sketched. The rodlike 
form and 10% extension (compared to dsDNA) of intermedi- 
ate DNA should be noted. 



DNA helices (Dff 6 = 3.2 • 1CT 11 m 2 /s) or ssDNA coils, 
and therefore can be representative only of intermediate 
DNA. Thus, the fraction of intermediate DNA im(c) as 
compared with fraction 1 — im{c) of native dsDNA across 
A and B regions is calculated in the following manner 



im(c) 



(D 



146 



(7) 



In Fig. [7] it is apparent that as DNA concentration is 
lowered dsDNA helices (open circles) give way to inter- 
mediate DNA (black circles). 



On lowering DNA concentration dsDNA disappears 
from population already in region B and only intermedi- 
ate DNA and separated strands of ssDNA could remain in 
solution. The variation of -D^gfc) across B,C,D regions 
may only be due to varying proportion of intermediate 
DNA and ssDNA coils, as they have different diffusion 
coefficients: = 2.9 • 1CT 11 and Df s 46 = 3.95 ■ 1CT 11 , 

respectively. Here DNA is in dilute regime and the diffu- 
sion coefficient of either conformation is not expected to 
depend on concentration. Thus, the fraction of interme- 
diate DNA im' '(c) as compared with fraction 1 — im'(c) 
of separated ssDNA strands across B,C,D regions is cal- 
culated in the following manner 



im'(c) 



(D 



146 



^146 > 



(8) 



In Fig. [7] it is apparent that as DNA concentration is 
lowered intermediate DNA (black squares) separates into 
ssDNA coils (open squares) , and in region D only ssDNA 
remains in solution. 

Finally, we note that im'(c) is complementary with 
im(c) data, and both indicate the same region c = 
0.05 — 0.3 mM where intermediate DNA is prevalent in 
the population. Interestingly, the fraction of open base- 
pairs in this region is about 50%, see FigEl This supports 
our initial proposition that we do not expect DNAs with 
varying proportions of open base-pairs. 

Half of the basepairs being open for intermediate DNA 
relates to the fact that DNA that we used has a quite ran- 
dom sequence and equal content of G-C and A-T pairs. 
Since A-T pairs have smaller binding energy [lj], it 
could be that only these break in the reported conditions 
and thus contribute to formation of intermediate DNA 
as a step in DNA melting process. 

Here, we remind that the counterion concentration is 
the only parameter that changes for this DNA and the 
reduction of screening due to decreasing concentration 
leads to DNA melting. Apparently, this parameter vari- 
ation is mild enough to allow for preferential breaking 
of A-T pairs, or, which is rather equivalent [13, 01) f° r ~ 
mation of bubbles in AT-rich portions of the sequence. 
That is, across our range of counterion concentrations 
1-0.01 mM the Debye screening length k^ 1 varies from 
about 10 up to about 100 nm. The corresponding ex- 
ponential factor exp(— nr) in the screened Debye-Hiickel 
potential at r = 2nm, characteristic distance of two re- 
pulsing phosphate charges would concomittantly change 
for about 20%. In other words, if the Coulombic re- 
pulsion at 1-2 mM counterion concentration was barely 
enough to break any basepairs, the reduction in screen- 
ing at lower concentrations would enhance the repulsion 
for 20% and start to break A-T pairs but not yet G-C 
pairs which are bound 50-100% stronger. 

Considering the above illustration we note that Chen 
and Prohofsky [5j calculate that below 1 mM added salt 
phosphate-phosphate force becomes insensitive to the 
salt concentration, and conclude that the effective poten- 



10 



tial has become a simple Coulomb potential. However, 
this is a long range potential and contributes to the en- 
ergy balance of the intermediate molecule specifically in 
the open, ssDNA sections, where the interaction occurs 
not only between complementary bases, but among at 
least several neighboring bases. That is, disrupted bond 
force constant is not zero. We emphasize that this ren- 
ders the calculation of free energy of such an intermedi- 
ate state nontrivial. Phosphate-phosphate contribution 
is non-additive, depends on DNA concentration, and on 
the size of the bubbles. As such it is not tractable and not 
predicted by present models for thermally induced DNA 
melting, where Coulomb potential was neglected, being 
screened (at most, Coulomb potential has only renormal- 
ized the hydrogen bonds). 

D. Concluding remarks 

In our scenario the intermediate DNA with about half 
basepairs open coexists with dsDNA helices or ssDNA 
coils in very low salt environment. We suggest that the 
low DNA concentration and very low salt conditions of 
this study allowed for DNA to open preferably along AT- 
rich portions of the sequence, without the strands sepa- 
rating further. These local lesions reflect on hydrody- 
namic properties of the whole molecule, i.e. its diffusion 
coefficient has been reduced measurably below the value 
for the original dsDNA helix. However, whether these 
DNA lesions remain stably open or it is a dynamical sit- 
uation where the bubbles open and close often enough to 
change the diffusion properties may not be judged from 
our result. Experiments capable of measuring the fluctu- 
ation rates at a given site within the sequence, as those 
mentioned in the Introduction (optical spectroscopy of G 
bases, FRET, FCS) should be performed on DNA frag- 
ments in the very low added salt conditions. Only in this 
manner the nature of intermediate DNA could be judged 
and then also the stability mechanism would be more eas- 
ily inferred. An interesting scenario would be the criti- 
cal slowing of fluctuations, a concept known from the 
studies of glass forming systems [46]. Finally, our find- 
ing of intermediate DNA provides another novel twist. 
That is, to explain the diffusion properties of intermedi- 



ate DNA in very low salt, the ssDNA loops within such 
a molecule have to be extended, and the molecule has to 
keep the rodlike form. This may be perceived contrary 
to the ususal notion that the fluctuational openings, for 
DNA in an added salt environment, remove the molecule 
away from the rigid rod form [2(], [2l|, [23J • 



V. SUMMARY 

In this work we have revealed the presence of novel 
intermediate DNA conformation at low DNA concentra- 
tions (c < 1 mM basepair) in very low salt (< 0.05mM) 
solutions at ambient temperature. The diffusion coefi- 
cient for DNA146 below 1 mM and above 0.05 mM (in 
DNA basepair) occurs to be below the values for both 
native dsDNA helices and ssDNA strands of equal poly- 
merization degree. Thus the population in the interme- 
diate regime is not a mixture of these two conformations. 
Intermediate DNA molecules, consisting of dsDNA seg- 
ments and ssDNA loops would have appropriate hydro- 
dynamical properties to explain for the diffusion data. 
The ssDNA loops are possibly due to DNA fluctuational 
openings (bubbles) that open and close at such rates, 
that they may be considered to have become stabilized 
since they influence the diffusion of the whole molecule. 
We suggest that the especially long-range Coulombic re- 
pulsion in the very low salt conditions of our experiment 
might allow for this effective stability. 
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